The linewidth of a non-Markovian atom laser 
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We present a fully quantum mechanical treatment of a 
single mode atom laser including pumping and output cou- 
pling. By ignoring atom-atom interactions, we have solved 
this model without making the Born-Markov approximation. 
We find substantially less gain narrowing than is predicted 
under that approximation. 
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An atom laser is a device which produces a coherent 
atomic de Broglie wave analogous to the coherent light 
wave produced by an optical laser. Unlike optical lasers, 
existing experimental atom lasers are not pumped KB. 
Consequently they do not exhibit gain narrowing, the 
phenomenon whereby the output linewidth is much nar- 
rower above threshold than below threshold. 

Atom laser models based on the Born-Markov approx- 
imations (BMA) predict gain narrowing, but they fail 
for a range of physically interesting parameter regimes 
In the BMA a laser, optical or atom, is described 
by a quantum optical master equation In the sim- 
plest case this models an oscillator subject to gain and 
loss, due to pumping and output coupling respectively. 
The linewidth is associated with the net dissipation in 
the system, which is the difference between the gain and 
loss. This difference decreases as the laser is pumped 
above threshold, giving rise to the gain-narrowed laser 
linewidth. This idealised Schawlow-Townes linewidth has 
its fundamental origin in spontaneous-emission-driven 
phase diffusion [JtJ. By analogy, Graham B has esti- 
mated the ultimate atom laser linewidth due to scatter- 
ing of thermally excited phonons. Previous studies based 
on the BMA have shown that certain pumping processes, 
such as evaporative cooling, severely broaden the atom 
laser linewidth beyond this fundamental limit [p|jL0|| . 

In the following we present a quantum mechanical 
analysis of the atom laser which does not make the BMA. 
We find that the gain narrowing is several orders of mag- 
nitude less than predicted using the BMA. Our methods 
may be useful for other non-Markovian systems, such as 
spontaneous emission in optical band-gap materials [ p"l| . 

Models based on the nonlinear Schrodinger equation 
are able to include interatomic interactions Jl2| — |l5|] . How- 
ever they implicitly assume that the lasing field is in a 
coherent state, which makes it impossible to calculate the 



linewidth of the resulting output, since no information 
remains about the quantum statistics. Interatomic inter- 
actions are difficult to include in a full quantum mechani- 
cal model, so our model assumes that they are negligible. 
This is an accurate description of very dilute systems. 

We have previously shown that a simple model of 
atom laser output coupling produces a nondispersing 
state which prevents a pumped atom laser from reaching 
a steady state ||. Adding gravity destroyed the non- 
dispersing state, overcoming the problem. Hence the 
atom laser modeled in this paper includes the effect of 
gravity on the output atom field. 

We model the atom laser by separating it into three 
parts. The lasing mode is an atomic cavity with large 
energy level separation. We assume that the cavity is 
single mode, with annihilation operator a and a Hamil- 
tonian H s . The external atomic field has a different elec- 
tronic state, so the atoms are no longer affected by the 
trapping potential. We model the external modes with 
the field operator ip(x) and the Hamiltonian H Q . The 
operators a and tp( x ) satisfy the normal boson commu- 
tation relations. The pump reservoir is coupled to the 
cavity by an irreversible process. At this stage, we will 
describe the pump by the Hamiltonian H p , which also 
couples the atoms from a pump reservoir into the system 
mode. The coupling between the lasing mode and the 
output modes is 
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where we have introduced the interaction picture opera- 
tors 



dx K(x)ipi(x, t), 
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J/h(H 3 +H p )(t-t ) 
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The shape of the coupling k(x) is determined largely by 
the spatial wavefunction of the laser mode. 

Using the unitary evolution operator corresponding to 
the interaction Hamiltonian Eq.(|l|), we find in the Heisen- 
berg picture, 

ip H {x,t) = ipi(x,t) - / ds F(x,t,s)a H {s), (2) 
J t 

where tpH(x,t) and au(t) are Heisenberg operators, and 
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F(x,t,s) = [ipi(x,t),tf{*)] 

■*«.(»)fefe<w}&,»)] 

dy K*(y)G(x,t,y,s) 
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where G(x, t, y, s) is the Green's function propagator due 
to the output Hamiltonian, H a , only. These functions can 
be written in closed form for several useful cases, includ- 
ing free space, free space with gravity, and a repulsive 
Gaussian potential. We may use Eq. (||) to calculate any 
observable of the output field, providing we know the 
complete history of the system, ajy(s). 

To calculate the output energy flux we transform our 
interaction Hamiltonian into the basis of the energy 
eigenstates of the output modes: H a = J dp fi w p c|cp, 
where c p is the annihilation operator associated with the 
eigenstate of H a that has a position space wavefunction 
u p (x) and energy fko p . Defining R,(p) — J dx u p (x) k(x), 
the output energy flux in terms of the two time correla- 
tion of the system is 



d{c\c 



dt 



2\R(p)\ 2 ^^ dse- iw * A V(i)a(s))Y (4) 



where Sft denotes the real part, and At = t — s. This 
assumes that at time t = 0, the output field was in the 
vacuum state. 

When the output field is in free space and the only 
term in H Q is the kinetic energy, then the eigenstates are 
the momentum eigenstates. In this case, R(p) is just the 
Fourier transform of k(x). When there is a gravitational 
field, the eigenstates are Airy functions with a displace- 
ment which depends on the energy: 



u p (x) — AfAi[f3(x — hujp/ (mg))] 



(5) 



where J\f is a normalisation constant, the length scale is 
given by (3 = (2m 2 g/h ) 1//3 , and m is the atomic mass. 
In this case ii(p) must be calculated numerically. 

Following Scully and Lamb we model pumping by the 
injection of a Poissonian sequence of excited atoms into 
the atom laser PJl^|. These atoms may spontaneously 
emit a photon and make a transition either into the atom 
lasing mode or into other modes of the lasing cavity. For 
simplicity, we consider a two-mode approximation. To 
obtain the pumping term, we consider the effect of a sin- 
gle atom injected into the atom laser, and then extend 
this to describe the effect of a distribution of atoms. This 
gives the master equation [ p"7| 



(P) 



pump 



rV[a)] 
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where r is the rate at which atoms are injected into the 
cavity, and n s the saturation number. The superopera- 
tors T> and A are defined by 



A[c], 



T>[c] — cpc^ 
A[c]p = - (c 1 cp + pc^ c). 



(7) 

(8) 



In our particular model n s depends on the ratio of the 
probability that an atom will spontaneously emit into the 
lasing mode to the probability that the atom will emit 
into another mode. 

If we make the Born approximation, trace over the out- 
put modes, and then make the Markov approximation, 
we can write the damping term of the master equation 

as a 



(^)damp = lBAiT>[a]a, 



(9) 



where jbm is related to the above threshold mean atom 
number n and saturation photon number n s by 



Ibm = r/(n + n s ). 
The full equation of motion is then 



P 



IBM 



dt V 2 J 

P« r/{2{n + n s ) + 1} 
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which has an error term proportional to r/n 3 if we as- 
sume that the system is close to a coherent state. The 
solution to this equation is an exponential decay, and the 
energy spectrum is therefore Lorentzian, with a width of 



Tbm = r/{m 2 ). 



(13) 



If we do not make the Born or Markov approximations, 
but we do assume that the trap population is localised 
around some (at this stage unknown) value n, we find 

| (at(t)) - (iu + P)(at( f )) + Jdx K*(x){4 >t ) (14) 

well above threshold. Note that we can no longer relate n 
directly to the physical parameters of the problem using 
Eq. (|l0|) , which has used the Born approximation. Since 
P determines n s , we seek an iterative method to produce 
a self-consistent solution. 

Under our assumptions the pumping is effectively lin- 
ear, so we may use the quantum regression theorem. Us- 
ing Eq.(^|), we may derive the following Volterra convo- 
lution type integro-diffcrcntial equation for the two time 
correlation function: 



d_ 

dr 



(at(t + r)a{t)) = {iuj + P)(at(* + T)a{t)) 

- / * +T du f*(t + t - u )(ot(«)a(t)), (15) 



where r > 0, and /(At) is the memory function 
f{At) = ( dx k{x)F{x,M). 



(16) 
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We have assumed that F(x, t, s) is a function of At = t — 
s. This again assumes that at time t = 0, the output field 
was in the vacuum state. This equation is not sufficient 
to specify the dynamics of the cavity, as it is only a single 
partial integro-differential equation in a two dimensional 
space. We also require the integro-differential equation 
for the intracavity number, which we can generate in a 
similar manner. Well above threshold, we obtain 



dt 



(ato>(t) 



ds 2${f*(At)(a^s)a(t))}. (17) 



These equations are difficult to solve in general, but 
can be solved in various limits. For example, if the kernel 
f{t) is a 8- function, as in the broadband limit of the 
optical laser, then the equations would become local, and 
the solution is an exponential. Although the broadband 
limit can be a good approximation for the atomic case as 
well [|, in general the atoms will disperse, which gives 
the system an irreducible memory. 

In order to produce analytical forms for the memory 
functions we assume that the coupling is Gaussian, and 
that there is no net momentum kick given to the atoms: 



K (x) = ^(2o|/7r) 1/4 exp[-(cr fc x) 2 ] 



(18) 



where Tiau is the momentum width of the coupling and 
7 is the strength of the coupling. In the presence of a 
gravitational field, V = mgx, the Green's function in 
Eq. (m can be found as a standard result |la,M : 



G{x,t,y,s) 



1 
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where A = h/(2m), and At = t — s. This leads to the 
following form for the memory function, 



/(At) 



7 exp( 



g 2 At 2 • 
' 32 A 2 a 2 ■ 



ex Pl 48 A 
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For coupling into free space, g = 0, and f(At) goes as 
l/VAt in the long time limit. There are no useful ap- 
proximations when g ^ 0. This is because the broadband 
limit of the integrals involving f{At) become unbounded 
in amplitude, and their convergence is due to their highly 
oscillatory nature. 

Eq. (jT|) and Eq. @ do not for m a standard pair of 
partial integro-differential equations. The derivative in 
Eq. (|l^) is only defined for t > 0, and so we cannot re- 
quire that the solution obey this equation through the 
whole domain of the integral. This means that Eq. @ 
cannot be integrated to find the solution, as we do not 
actually know the derivative at any point. Since the two- 
time correlation is Hermitian, we can we rewrite the in- 
tegral so that the domain remains in the r > plane, 



but we still do not have a continuously defined derivative 
along the length of the integral. 

We proceed by making an ansatz which uses the solu- 
tion of Eq. (|lj) which has been extended into the region 
r < 0. We then use the r > portion of this solution to 
substitute into the two time correlation in Eq. (|l7|). We 
introduce the function J it), which is the solution of the 
equation 



dJ(t) 
dt 



(iUJ + P)J(t) 



ds f*(At)J{s) 



This means that 



(a\t + T)a(t)) = {a!a)(t)J(t + r)/J(t) 



(21) 



(22) 



is a solution of Eq. (|15|) with the correct initial condition 
at r = 0. We then substitute this result into Eq.(p"7|): 



dt 



: . M( . )B (3«*jp). (23) 



Solving these two equations gives the two time corre- 
lation for the lasing mode, from which we may find the 
properties of the output field. It is only consistent with 
our linearisation of the pumping if the number of atoms 
in the trap, n, converges to the value which originally 
produced the parameter P. Since we require P to gen- 
erate the solution, and n is simply the long time limit of 
Eq. (p3|), the effective free parameter is n s . This thresh- 
old parameter must be much smaller than n, so we search 
for a value of P which gives the result n 3> n s . 

Once it is established that Eq. (23) is approaching a 
stable steady state, a fast way of finding it is to set the 
derivative to zero, and assume that (a'a)(s) = n over the 
support of the kernel. This gives 



V J ( s ) J 



(24) 



We transform Eq. ( pi| ) into a rotating frame by intro- 
ducing the function N(t) = J(t) exp(— iu) t — Pt) and 
rewriting it in terms of N(t): 



dN(t) 
dt 



ds H(At)N{s) 



(25) 



where H(At) = f*(At) exp(-[iw + P]At). An effec- 
tive numerical method for solving this equation is direct 
integration using a second-order algorithm for both the 
integration and the calculation of the integral to find the 
derivative at each timestep [M . 

We now calculate the output properties of an atom 
laser. We use the trap frequency uo — 2ir x 123 Hz 
0, an atomic mass of m = 5 x 10~ 26 kg, a gravitational 
acceleration of g = 9.8sin(0.18) m s~ 2 , and the coupling 
given by Eq. ( p"8| ) with a momentum width — 4.4 x 10 5 
m _1 . We use a damping constant of 7 = 2.0 x 10 4 s -2 , 
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and a threshold of n s — 47. As the pumping rate r 
increases, the steady state number of atoms in the cavity 
n increases, and the modulus of the two time correlation 
decays more slowly. The energy spectrum of the output 
flux is proportional to the Fourier transform of the two 
time correlation through Eq. (|4|), so the laser shows gain 
narrowing. 

In Table 1 we show the results of these calculations. 
The linewidth of the output flux, T, is calculated di- 
rectly from the two time correlation using Eq. (Q). This 
is compared to the linewidth given by the Born-Markov 
approximation Tbm, Eq. (|l3|). 

The atom laser linewidth predicted by our theory is 
several orders of magnitude larger than that predicted 
using the BMA result. We plot the spectral flux corre- 
sponding to three of these pumping rates in Fig. [ij The 
vertical scale is normalised to the peak height for each 
plot so that the width of the spectra can be easily com- 
pared. The spectra are almost Lorenztian, but have a 
drastically different width and are slightly shifted com- 
pared to the results under the BMA. 

Our model does not include atom-atom interactions, 
and therefore only works when the atomic field is very 
dilute. The advantage of our model is that, unlike mean 
field models, we have avoided making the approximation 
that the lasing mode is perfectly coherent. Future work 
will involve generalising this atom laser model to include 
atom-atom interactions. 
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TABLE I. Linewidths and cavity atom number as a func- 
tion of pumping rate r. 
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1.1 
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0.035 


0.00062 
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FIG. 1. The spectra of the output energy flux for three 
different pumping rates. The dotted, dashed and solid lines 
represent r — 40s -1 , r — 80s -1 and r = 800s -1 respectively. 
Other parameters are given in the text. 
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